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ABSTRACT 

Several works have reported changes of the Sun's subsurface stratification 
inferred from /-mode or p-mode observations. Recently a non-homologous varia- 
tion of the subsurface layers with depth and time has been deduced from /-modes. 
Progress on this important transition zone between the solar interior and the ex- 
ternal part supposes a good understanding of the interplay between the different 
processes which contribute to this variation. This paper is the first of a series 
where we aim to study these layers from the theoretical point of view. For this 
first paper, we use solar models obtained with the CESAM code, in its classi- 
cal form, and analyze the properties of the computed theoretical /-modes. We 
examine how a pure variation in the calibrated radius influences the subsurface 
structure and we show also the impact of an additional change of composition on 
the same layers. Then we use an inversion procedure to quantify the correspond- 
ing /-mode variation and their capacity to infer the radius variation. We deduce 
an estimate of the amplitude of the 11-year cyclic photospheric radius variation. 

Subject headings: Sun: helioseismology — Sun: oscillations — Sun: activity — 
Sun: interior — Sun: model — Sun: radius — Sun: luminosity — Sun: variability 



1. Introduction 

Helioseismology has been extremely useful to probe the internal structure of the Sun 
but two crucial regions need to be improved for a good understanding of the time evolution 
of the solar activity: (1) the solar core for a proper description of the transport of momen- 
tum implied by rotation, gravity waves and magnetic field along the evolution and (2) the 
subsurface layers which correspond to the transition region between large and small scale 
dynamics evolution. 
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This transition zone is important to study because it couples the internal dynamics to 
the dynamics of the solar atmosp here and it plays a crucial r ole in the emergence of the space 
weather science (see for example iRozelot &: Lefebvrd (l2006l )). Figure [T] sketches a schematic 
view of these subsurface layers above 0.96 Rq, which are called the leptocline region (f r om th e 
greek "leptos"= thin, "klino"= tilt). This term was proposed by iGodier fc RozelotI (120011 ) . 
who computed the solar oblateness and showed a curvature change in this zone which was 
interpreted as the presence of a double layer. Just below the surface, there is a change in 
radial and latitudinal rotation (IBasu. Antia fc Tripathyl 1 19991 : ICorbard fc Thompson! |2002| ) 
and the treatment of the superadiabiatic region supposes a proper description of the con- 
vection (presently described by the mixing length parameter) and of detailed molecular and 
atomic opacity calculations. One needs also to add the turbulent pressure and the emer- 
gence of local magnetic field to these processes, then below, hydrogen and helium pass from 
neutral to partially ionized and then totally ionized in a region where the magnetic pressure 
cannot be igr iored. The mean and varying magnetic field amplitudes have been tentatively 
extracted by iNghiem et al.l (120061 ) from the analysis of the low degree acoustic modes. The 
proper interaction between these different physical processes is not yet included in stellar 
evolution models and could lead to some variation of the solar radius along the solar cycle. 

The observation of these subsurfac e layers is not so easy. The difficulty com es from 
the fact that the high-degree p-modes ( iKorzennik. Rabello-Soares fc Schou I |200J) have a 
small lifetime and are largely perturbed by the turbulent motions and the emerging mag- 
netic field. But impo rtant progress has appeared recently than ks to the analyses of long 
series of seismic data. iRabello-Soares. Korzennik &: Schoij (120081 ) estimate the variations of 
the high-degree p-modes fre quencies over the solar cycle and using ring diagram analysis 
Basu. Antia fc BogartI (120071 ) extract now latitudinal and temporal variations of the sound 
speed or the density along the Hale solar cycle (at least its last half cycle). 
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In parallel, iLefebvre fc Kosovichevl (120051 ) and iLefebvre. Kosovichev &: RozelotI (120071 ) 
reported changes with the solar cycle of the solar subsurface stratification, inferred from 
inversion of SOHO/MDI /-mode frequencies. They notice a different behavior for the layers 
around 0.99 Rq. Between 0.97 and 0.99 Rq, it seems that the position of the layers varies 
in phase with the solar cycle, whereas it appears opposite in the upper part, above 0.99 
Rq, where the variability is in antiphase. So, it seems that the most external layers of the 
Sun shrink during the ascending phase of the solar cycle and relax after the maximum, but 
these changes are not uniform with depth. At the surface, this result is coherent with the 
observation of the photospheri c radius variatioris. Th e observed series of ground-based mea- 



surements (iLaclare et al. 



1996 



Reis Neto et al.ll2003l ) and the measurements aboard strato- 



spheric balloons ( Sofia. Heaps &: Twigg 1994j ) suggest a reduction of the solar photospheric 



radius at the maximum of the cycle. The balloon variation estimates are much smaller than 
the ground-based observations , pollu ted by the atmospheric turbulence. They are consistent 
with results from iKuhn et al.l (120041 ) who reported no evidence of solar-cycle visible radius 
va riations between 1 996 and 2004 larger than 7 mas. But they disagree with computations 
of ISofia et al.l (120051 ). These authors predict also non- homologous subsurface stratification 
changes with the 11 year solar cycle but with an amplification up to a factor 1000 from the 
depth at 5 Mm to the surface. It leads to a variation of the solar radius up to 600 km along 
the cycle as shown in Fig. 2 of iLefebvre. Kosovichev fc RozelotI ( 120071 ). Such a variation 
seems extremely large. 

So the study of the detected /-modes is v ery useful in the present context even the 
interpretation of /-modes has appe ared puzzling (jPziembowski. Goode fc Schoij|200ll : lAntia 
20031 : iDziembowski fc Goodell2004l ) due to some apparent lack of sensitivity of these modes 
at the real surface. We develop here a theoretical approach to vahdate the /-modes inversion 
procedure used previously and applied here to some specific cases that we know perfectly. 
This paper is the first of a series where we estimate the impact of a change of radius and 
composition on the subsurface layers and on the /-mode frequencies using classical solar 
models including a detailed microscopic description of these layers. In the present work, we 
do not try to justify the origin of the variation of the radius. The next step will be to develop 
more complex models including magnetic field and differential rotation. Our general aim is 
to properly qualify the variabilities of /-modes and solar radius that we will study with 
the SDO and PICARD missions (Kosovichev et al. 2007 and Thuillier, Dewitte & Schmutz 
2006, respectively). We will also investigate the capability of the /-modes to estimate the 
solar photospheric radius variability. This study will contribute also to clarify the notion of 
"solar radius" to bridge two different c ommunities and obtain without ambiguity a unique 
definition of this fundamental quantity (IHaberreiter. Schmutz fc Kosovichevl l2008l ). 



After a description of our study (Section [2]), we describe in Section [3|, the models used 
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and examine the subsurface layer changes on different variables produced by a change in 
the solar radius and composition and we calculate the corresponding /-mode frequencies. 
Section H] is devoted to the use of these frequencies to infer the change in the position of the 
subsurface layers. The validation of the procedure, an estimate of the present solar cyclic 
radius variation and the perspectives of this work will conclude the paper in Section [51 



2. A new estimate of the solar cycle radius change using models and f-modes 

predictions 

The present the oretical work focus on the solar layers located above 0.96 R0. The 



previous approach of iLi. Basu fc Sofia! (120031 ) tried to interpret the behavior of the seismic 
acoustic observations by introducing some magnetic effect in solar models to predict the 
radius variation over the cycle. In the present work, we concentrate on the /-mode frequen- 
cies. We use known models and perform inversions that we can verify and compare with 
the present observational /-mode frequency variations. For this first paper, we stay in the 
classical approach and only focus on basic quantities which may vary along the solar cycle 
without introducing the physical process at the origin of the corresponding variations. 

The next obiectiy es will be to introduce the different magnetohydrodynamical actors 



f lMathis fc Zahn 



2005[). Turbuleii ce. rotation and the local effect of magnetic field with a 



decent topology (IDuez et al.ll2007l ) must be introduced together in a more sophisticated way 
than it is generally done. But we would like first to study the potential of the /-modes in 
some well known cases and then show that it is possible to deduce an order of magnitude of 
the solar radius variation from the present variation of the observed /-mode frequencies. 

The general idea is to progress step by step on the sensitivity of the /-modes to the 

> using the same formalism than the one used in 
present p aper we do not yet impose any k ind o f 



Lefebvre & Kosovichev 


( 


2005). In the 


magnetic field like in 


Li. Basu & Sofia 



but we estimate the local changes to decouple the origin of the different effects. We examine 
the theoretical /-mode frequency changes coming from known solar models which mimic a 
simple expansion or a change in composition. We deduce, for the first time from /-modes, 
an estimate on the quality of extraction of the photospheric radius change and an order of 
magnitude of its Hale cycle change. 
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3. Model computations 



We use the CESAM evolution code (IMorell 119971 ) to calculate several models of chosen 



solar radius R, luminosity L and surface abundance (Z/X). The different computed cases are 



organized in 3 sets, and hsted in Table [H Our reference m odel Mi ( iTurck-Chieze. Couvidat fc Kosovichev 



200ll : ICouvidat. Turck-Chieze. S.. fc KosovichevI |2003|) with Rq = 6.9599 x lO^km, Lq 



3.8460 X lO^^erg/s and Z/Xq = 0.02447 is the seismic solar model (SeSM) built to repro- 
duce the observed sound speed profile. 

Set 1 is composed of models M2 and M3 converged with a radius fixed respectively 
at = 1.0002 and = 0.9998 and the composition fixed at the value of the reference 
model, but where the mixing-length parameter a is a free parameter (the luminosity is not 
constrained in models M2 and M3). The radius difference of 2 x 10^^ corresponds to a 
difference of about 140 km. 

Set 2 is composed of six models Mj, for i = 4. ..9 built with the two free parameters Yi 
the initial composition in helium and a; in these models we impose a change in the radius 
and luminosity, and the same ratio Z/ Xsuri of the final superficial chemical composition than 
that of model Mi. 

Finally set 3 consists of five models Mj, for i = 10. ..14 which differ from set 2 by 
the value of Z/Xg^^^, and includes a model where R = Rq and L = Lq (M14). A change 
of composition is interesting to study because the presence of varying magnetic field in 
the subsurface layers may lead to an apparent change of hydrogen/helium thermodynamic 
conditions in solar models with some impact on the density and pressure of the subsurface 
layers. 

The goal was to mimic variations of the radius and luminosity which could be extrapo- 
lated to variations during the 11-year solar activity cycle. For instance, M4 and Mg have a 
luminosity varying in the same way than the radius, whereas Mq and M7 present luminosity 
and radius varying in an opposite way. The relative variations 2 x 10""^ for the radius are 
chosen larger than the supposed observed ones in order to avoid numerical accuracy prob- 
lems and to better see their influence on the subsurface layers. In the different cases, the 
models converge to get the calibrated radius within a precision of 10^^ that means within 7 
km. The relative change in the luminosity of 10~'^ is similar to the observed solar irradiance 
variation during the 11-year cycle. The way we introduce it in solar models leads generally 
to a variation of luminosity on the production of energy instead on the superficial layers. 
So we see only an indirect effect of the luminosity change on the subsurface layers through 
changes in temperature and radius, their effect is not negligible for models 10 and 13. 
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3.1. Changes of the subsurface layers 

We estimate how the subsurface layers react to the induced changes. For that, we 
estimate the differences of the following variables between the different models and the 
=ei=mic modeQ the ma=s ,n, the temperature T, the density p, the pressure p, the =ound 

speed c, the adiabatic exponent Fi = (^||^)adia5 the density scale height Hp = — (^^^) , 

the pressure scale height Hp = — (^^) ^, the radiative temperature gradient Vrad = 

the real temperature gradient Vreai = iiiin(Vrad) Vadia), the Rosseland opacity coefficient k, 

and the gravitational energy Eg. 



^Note that each fractional radius Xi is reported to the reference radius, i.e. cc"'^™ = Xi Ri/Rref to get all 
models refered to the reference model. 
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Table 1: List of the calculated models. 



Model 








Yi 




a 


R/Rq 






TefAK] 












Seismic model 














1 


0, 


.70642 


0. 


,27468 


2, 


.03990 


1.0000 


1 


.0000 


0, 


.02447 


5777. 


,52 














Set 1 














2 


0, 


.70642 


0. 


,27468 


2, 


.03751 


1.0002 


1 


.0000 


0, 


.02447 


5776, 


,97 


3 


0, 


.70642 


0. 


,27468 


2, 


.04260 


0.9998 


1 


.0001 


0, 


.02447 


5778. 


,30 














Set 2 














4 


0, 


.70632 


0. 


,27478 


2, 


.04033 


1.0002 


1, 


.0010 


0, 


.02447 


5778. 


,40 


5 


0, 


.70642 


0. 


,27468 


2, 


.03757 


1.0002 


1 


.0000 


0, 


.02447 


5776. 


,98 


6 


0, 


.70651 


0. 


,27458 


2, 


.03463 


1.0002 





.9990 


0, 


.02447 


5775. 


,52 


7 


0, 


.70634 


0. 


,27476 


2, 


.04508 


0.9998 


1 


.0010 


0, 


.02447 


5779. 


,56 


8 


0, 


.70644 


0. 


,27466 


2, 


.04214 


0.9998 


1 


.0000 


0, 


.02447 


5778, 


,10 


9 


0, 


.70653 


0. 


,27457 


2, 


.03937 


0.9998 





.9990 


0, 


.02447 


5776. 


,67 














Set 3 














10 


0, 


.69772 


0. 


,28186 


2, 


.06148 


1.0002 


1 


.0010 


0, 


.02676 


5778, 


,42 


11 


0, 


.69791 


0. 


,28166 


2, 


.05575 


1.0002 





.9990 


0, 


.02676 


5775. 


,53 


12 


0, 


.69774 


0. 


,28184 


2, 


.06624 


0.9998 


1 


.0010 


0, 


.02676 


5779. 


,57 


13 


0, 


.69793 


0. 


,28164 


2, 


.06048 


0.9998 





.9990 


0, 


.02676 


5776. 


,69 


11 


0, 


.69783 


0, 


.28175 


2 


,06098 


1.0000 


1 


.0000 


0, 


,02676 


5777, 


.00 



Note. Model 1 is the reference model (Turck-Chiczc et al. 2001; Couvidat, Turck-Chieze and Koso- 
vichev, 2003), Xj and Yi are respectively the initial composition in hydrogen and helium, a is the mixing- 
length parameter, R and L are the radius and luminosity, Z/X^^^f is the final superficial abundance and 
Te// the eS'ective temperature. 



- 9 - 




Fig. 2. — Difference for the physical variables between models of set 2 calibrated with 
different combinations of radius and luminosity (see text) and the seismic solar model. The 
curves of set 1 are exactly superimposed with those of set 2. M2 (or M3) and M5 (or Mg) 
looks hke M4 and Mg (or M7 and Mg). 
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Fig. 3. — Difference for the physical variables between models of set 3 and the seismic solar 
model (same principle than in Figure [2] but including an additional effect in the composition). 
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Figures [2] and [3] illustrate the changes with depth of these variables in comparison with 
the seismic one. We focus our study to the zone located above 0.96 Rc:^, that is the zone where 



changes in the subsurface stratification have been found by iLefebvre &: KosovichevI (120051 ) 
along the solar cycle. The upper limit is fixed at 0.998 Rq, beyond which the superadiabatic 
zone extends, the turbulence acts strongly and the rotation changes quickly. We have chosen 
the zone where the /-modes have a good sensitivity. At a first glance, we notice that most 
of these variables present non negligible variations in the studied region. 

We first comment on the common trends of both figures: 

• The variations shown on the different panels are dominated by the change in radius 
imposed in the calibration and at the second order by the change of composition for 
set 3. 

• The way we introduce the change in luminosity acts on the nuclear burning layers and 
has a negligible effect on the subsurface layers. So for clarity, we have not plotted the 
models 2 and 3 of set 1, nor 5 and 8 of set 2, in fact their curves are aligned with 
the models having the same radius in Figure O In contrast the observed luminosity 
variation of 10~^ along the solar cycle is most likely coming from the very external 
layers. In fact, the variation of radius and of the photospheric temperature of our 
models produces a change of luminosity determined by the Stefan's law. It is generally 
too small to exhibit any structural effect, except for models 10 and 13 where the 
variation of luminosity coming from the variation of radius and temperature is near 
from the imposed variation of luminosity. 

• The behaviors of c. Hp and Hp are very similar i.e. a bump around 0.99 Rq with 
opposite variations between models of different radius; this bump exists also for the 
temperature. This position corres ponds to the transition between neutral He and He"*" 



as discussed in iLopes et al.l (119971 ) 



• The differences of Fi and Vreai have similar variations and present a double peak near 
0.99 -R0 with almost equal amplitude. In addition to the bump at 0.99 Rq, the second 
bump could be due to the transition H+/H. It is reasonable that in this convective 
region, the gradient of the structure follows the adiabatic exponent. 

• The variations of k, and Vrad present a sign change at 0.99 Rq, connected to the vari- 
ation of the opacities in the region where the light elements are partially ionized. The 
variation of the pressure and density modifies the corresponding opacity coefficients. 



Figure [3] differs nevertheless from Figure [2j Figure [2] shows two symmetrical groups of 
models depending on the value of the calibrated radius. In the set 3 another change comes 
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from the modification of the heavy elements contribution (about 10%) which induces a small 
change in helium (about 3%) and hydrogen. In order to separate properly the effect of radius 
from the effect of composition, we have calculated model (M14) that includes only the change 
in composition. We show in Figure [3] the changes induced by the composition effect. As 
an evident consequence, we loose the symmetry between the group of models with a greater 
radius compared to the models with a smaller one, this is mainly visible on the temperature 
and Fi differences. 

These first behaviors indicate that the subsurface layers above 0.96 Rq are significantly 
affected by a change in radius and composition. To go further, we shall estimate the differ- 
ences in /-mode frequencies issued from these models. 



3.2. Calculation of the theoretical /-modes for the different models 

We compu te for each model the theor etical /-mode frequencies using the oscillation 



code ADIPLS (iChristensen-Dalsgaardl Il982l ) to see the impact of the different changes on 
these quantities. Figure H] shows the relative difference of frequencies between each model 
and the reference model Mi versus the corresponding absolute frequencies. 

We note first that the /-mode frequencies of a model with a larger radius are smaller 
than those of the reference model. The change of 2 x 10~^ in radius leads to a relative change 
of about 3 X 10"'' on /-mode frequencies. 

The left panel of Figure 4 corresponds to sets 1 and 2. The difference associated to 
M4 and Mg is slightly bigger in absolute value from the difference associated to M7 and 
Mg due to the value of the reached real radius: +135 km for M4 and Mg and -143 km 
for M7 and Mg (this difference is within the precision imposed of the model radius). The 
frequency dependence has almost a fiat behavior, whereas this is not the case for set 3 in 
the right panel where the relative difference has an extremum around ~ 1300 fiHz due to 
the additional effect of the change in composition, which produces non-uniform changes as 
already discussed. 
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Fig. 4. — Relative difference of /-modes frequencies between models of set 2 (left panel) or 
models of set 3 (right panel) and Mi versus the absolute frequencies. We do not plot the 
difference between the two reference models M14 and Mi because this difference is of the 
order of 10~^. 
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4. Numerical inversion of theoretical /-mode frequency to infer subsurface 

stratification changes 

In this section we do the inversion of the /-mode frequencies for degrees / between 100 
and 300, in a range slightly b igger than the one chosen pr eviously with the corresponding 
solar observational quantities f Lefebvre fc Kosovichevl[2"o05l ) to see the reproductivity of the 
radial variation for known models, to qualify the procedure and to extrapolate a new estimate 
of the radius variation along the solar cycle. 



We use the same formalism as in iLefebvre fc KosovichevI (120051 ) to infer the changes in 
the position of subsurface layers from the /-mode frequency variations. A relation between 
the relative frequency variations 5v/v for /-modes and the associate d Lagrangian perturba- 
tion o f the radius 5r/r of the subsurface layers has been established by lDziembowski fc Goode 



31 



where I is the degree of the /-modes, / is the moment of inertia as classically defined in 
Dziembowski fc Goodd (120041 ) . u is the angular frequency of the eigenpulsation (cu = 27rz/) 
and g is the gravity acceleration. The validity of this equation is lir nited to the case where 
any m agnetic field effect is explicitly introduced in the equations (IDziembowski fc Goode 



20041 ) which is the case in this work. It also supposes that, if y and z are respectively the 
vertical and horizontal eigenfunctions, the property of /-modes y'^ = l{l + Vjz^ is satisfied. 

This equation allows us to obtain 5r from 5v/ v. For these inversions, we used Mi as the 
reference model and a standar d Regularized Least-Square technique, since Equation [1] defines 
an ill-posed inverse problem (ITikhonov fc Arseninlll977l )). The validity of this equation for 
the different cases analyzed here will be discussed in the last section. 
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x=r/R degree I degree I 



Fig. 5. — Left panel: Radial variation Ar as a function of the fractional radius x = r/Rref, 
obtained by inversion of Equation [T] with frequency variations represented in the left panel 
of Figure m for set 2. The error bars are the standard deviation after averaging over a set of 
random noise added to the relative frequencies reconstructed in the right panel. Right panel: 
For each model, relative variation of frequencies Au/uref as a function of the degree /. The 
points with error bars (that have been arbitrarily added to the data, such as cr = 0.5 x 10~^) 
represent the input data, and the solid curves are the result of direct integration of Eq. [H 
providing the solutions in the left panel. 
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Fig. 6. — Idem than Figure [5] for models of set 3. Error bars have been arbitrarily added to 
the data, such as cr = 1 x 10~^. Note here the behavior of the curves near x = 0.99 reflecting 
a non-uniform change in the subsurface stratification. For clarity we do not plot the relative 
variation of frequencies for model M14 that is 10 times smaller than for the other models. 
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Figures [5] and [6] show results from inverting Equation [T] for each computed model. These 
figures also show the frequencies computed by integrating the solutions using Equation [H 
and how well they match the model frequencies within the error bars. The uncertainties were 
set arbitrarily to 0.5 x 10~^ and 1 x 10~^ respectively for the need of the inversion process 
and the will to mimic real variations. This fact guarantees the good quality of the inversion. 
The main characteristics of our solutions are: 



The inversion solutions plotted in Figure for set 2 lead to an almost uniform radial 
variation. The amplitude of the radial variation is similar to the one imposed at the 
surface, as expected from the relative variations of frequencies in the left panel of 
Figure HJ i.e. about ±140 km with a sign respecting the nominal one. Nevertheless, 
there is a small difference near the surface (see also next section), which is due to the 
poor spa tial resolution of the ke r nels a t the surface. This problem has already been 



raised in iLefebvre fc KosovichevI (120051 ) where the shape of the kernels is shown. 



• The solutions for set 3 shown in Figure E] are slightly different. Below O.98i?0, the 
variation is constant but closer to the surface there are non-monotonic changes in the 
stratification with a bump centered at O.99i?0. The uncertainty of the localization of 
this bump is governed by the characteristic width of our kernels that is about 0.005Rq. 



These variations are different in shape and in amplitude to those found by lLefebvre fc Kosovichev 



(120051 ). This is not surprising because we do not yet introduce the dynamical processes 



which generate the solar cycle. 

Since the only difference between some models of set 2 and some others of set 3 is the 
different value of {Z/X)^j, this bump is clearly a consequence of the introduction of 
this change in composition affecting the subsurface layers through the pressure and 
mainly the equation of state. We note that the value of AR is not any longer about 
140 km but about 110 km (135-20 km) for models Mio and Mn and -160 km (142+20 
km) for Mi2 and M13 at O.96i?0, this is easily explained by looking at the AR between 
the two reference models of about 20 km at this depth. It is clear that a change of 
composition has an effect largely below the superficial layers. 



In fact the inversion process implies that the solution is not unique. By adjusting the 
regular isat ion parameters, we can find another profile with two bumps for Ar with also a 
very good fit to frequency. That solution has not been retained because (i) the variations at 
the surface are not in agreement with the input radius variation as cited in Table [T], (ii) the 
error bars are bigger and (iii) the curves are more oscillating. Moreover, we cannot presently 
give a physical interpretation to this solution. On the contrary, for the first solution, we 
discuss in section |5] how we understand the present inversions. 
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Discussion 



Equation 8 of iDziembowski fc Goodd (120041 ) supposes that the radial position r used in 
Equation [1] is the lagrangian radius, /-modes are surface oscillation trapped waves, so each 
mode oscillates around an equilibrium radius. We relate in this section this radial position 
to that of the structure. 



5.1. Validity of the procedure of inversion 



The extraction of the solar subsurface structure and the photospheric radius variation 
along the solar cycle is not so easy to determine and several papers have been published wit h 



different conclusions (IDziembowski. Goode fc Schoul l200ll : lAntial l2003l : ISofia et al.l |2005| ) . 



The most recent work of IDziembowski &: Goodd (120041 ) shows the difficulty to obtain a 
general expression for the inversion, so this work allows to test such a procedure. We verify 
in this section if the displacement Ar of Equation [1] corresponds to the radial displacement 
at a given mass for the different cases studied. 
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Fig. 7. — Left panel: Radial displacement Ar of a given m(r), p(r), c,Hp, Ti and k defined 
in the reference model (Mi) for models of set 2 (the curves of set 1 models superpose the 
previous ones). These curves have to be compared to the left panel of Figure [H Right panel: 
Idem with models of set 3. 
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Fig. 8. — Radial displacement Ar issued from the inversion (same curves than in the left 
panel of Figure [6] added with the associated displacement at Hp constant as plotted in the 
right panel of Figure [7] (here dashed curves associated at the same color for the corresponding 
model). Note the good superposition of the different curves. 
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With this purpose, we calculate for each set (sets 1 and 2 lead to the same result) the 
radial displacement of a layer corresponding to a given value of the structural quantities, i.e. 
Ar = (r^r — ri)gj^=gj where q represents a model structure quantity like m, c. Hp or Fi, and 
is the model number. 

Figure [7] shows the panels representing Ar at constant m, c. Hp and Fi. For sets 1 or 
2, we have the same behaviour for all the quantities and Equation [1] can be applied without 
any doubt. For the set 3, the radial difference Ar, computed by the procedure described 
above, differs from one quantity to another. The only curve similar in amplitude and in 
shape with the one computed by inversion of /-mode frequencies in Figure [6] is the curve 
computed at fixed pressure height Hp. Effectively Figure [8] shows an excellent superposition 
of the corresponding curve and the radial displacement Ar issued from the inversion. This 
agreement, in shape and amplitude, demonstrates that the validity of equation 1 in this case 
requires to associate the radial displacement to this quantity. It is in fact not surprising 
because the /-mode frequencies are naturally sensitive to the pressure scale height. 

We note also that the inversion is not able to reproduce precisely the behavior very near 
the surface due to the lack of sensitivity of the /-modes. But we can deduce from this study 
an estimate of the uncertainty on the determination of the photospheric radius variation: we 
estimate it of the order of 15%. 



5.2. A new prediction on the solar radius variation along the solar cycle 

We show in this paper that a pure variation in the solar radius produces changes in the 
subsurface layers that a r e in t he same zone than that studied with the observed /-modes by 



Lefebvre fc KosovichevI (120051 ). Such a change is characterized by changes in the subsurface 
stratification and more exactly variations in the computed position of the layers. We show 
that the changes detected by /-modes are physically related to the variation of the pres- 
sure. In our modelling we are able to produce non-uniform changes by slightly modifying 
the chemical composition below the solar surface. This is not excluded if one introduces 
dynamical processes in keeping the same constraints on luminosity and radius. The two 
studies show different stratifications of the outlayers but with similar effects on the /-mode 
frequencies. Considering larger effects on the /-mode frequencies than in the solar case, we 
propose to deduce from this study an estimate of the variation of the solar radius along the 
cycle. From Au/u of about 3 x 10^^ for a Ar of about 140 km, we lead to a solar radius vari- 
ation along the solar cycle of about 7± 1 km for the observed variatio ns in frequency of about 



1.5 X 10 . This result is consistent with the radius extrapolated by lLefebvre &: Kosovichev 



(120051 ) (see respectively Figures 1 and 3 of this article). It is interesting to note that this 
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radius variation estimate is also compatible (in order of magnitude) with the observed low 
degree ac oustic mode varia t ion of about 0.4 uHz alorig the solar cycle of the low degree 
p-modes (IFossat et al.l 119871 : IChaplin. Elsworth fc IsaakI 120011 ). 



5.3. Perspectives 



This study represents the first step on the way to explain the variations of the subsurface 
stratification along the solar cycle and reinforces the interest of the /-modes. The next step 
will be the development of models including magnetic field and other dynamical processes, 
which will permit a more realistic study. Nevertheless it will suppose also an improved 
expression for the inversion of /-modes which is not yet obtained. We believe that the 
introduction of a magnetic f ield w ill influence the stratification of the subsurface layers. 
The study by iNghiem et al.l (120061 ) pointed out non-uniform radial changes, depending on 
the importance of the magnetic pressure. Moreover, it will be also interesting to take into 
account the differential rotation and consequently some asphericity to confront them to the 
subsurface latitudinal stratification over the 11-year cycle. 

We would like to emphasize that a better knowledge of these subsurface layers will 
contribute to our understanding (i) of the dynamics of the 11-year solar cycle and (ii) 
of the Sun-Earth relationship for space climate. This supposes in parallel the simultane- 
ous measurement of the radius and frequencies change. This study is in the framework of 



coming future space missions: SDO (see http://sdo.p:sfc.nasa.g:ov/)) a nd PICARD (see 



http : //smsc . cnes . f r/PICARD/ and iThuillier. Dewitte fc Schmutzl (120061 )). 



The DynaMICCS/HIRISE perspective (ITurck-Chieze et al.ll2006l . l2008[ ). proposed in the 
framework of the ESA Cosmic Vision 2015-2025, could bring even more data on all sources 
of internal variability by putting together instruments which will follow seismically all the 
layers down to the core together with instruments measuring the variability of the above 
atmosphere. 



We would like to thank two anonymous referees who have helped us to clarify our results 
and to enrich this paper. S. Lefebvre was supported by a CNES/GOLF postdoctoral grant 
in SAp/CEA. 
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